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ABSTRACT 

Recently, OGLE team have reported clear quasar microlensing signal 
in Q2237+0305. We have analyzed the microlens event of "image C" by 
using their finely and densely sampled lightcurves. From lightcurve fitting, 
we can unambiguously set the source size of 0.98 Einstein Ring radius 
as a conservative limit. This limit corresponds to 2000 (AU), if we adopt 
M\ ens ~ O.1M obtained by a recent statistical study of mean mass of lens 
object. This gives a clear evidence for the existence of an accretion disk in the 
central region of the quasar. 



Subject headings: accretion, accretion disks — galaxies: active — gravitational 
lensing — quasars: individuals (Q2237+0305) 



1. INTRODUCTION 

It is widely believed that central origins driving activity of quasars and AGNs (Active 
Galactic Nuclei) is an accretion disk surrounding a 1O 6 ~ 9 M SMBH (supermassive black 
hole). To find a direct evidence for this general belief is one of the most exciting subjects 
in the current research in astronomy and astrophysics, but unfortunately, the expected 
angular size of accretion disks is too small ( ^ 1 fias) to directly resolve spatially by using 
present observational instruments. For this reason, a proof of SMBH hypothesis remains as 
an unsolved problem in fields of quasar/AGN. Such a situation will not alter in near future. 
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However, there is a strong tool to make it possible. That is so-called "quasar 
microlensing" . Following the first report of a detection of quasar microlensing in 
Q2237+0305 ("Einstein Cross" or "Huchra's lens") by Irwin et al. (1989) and subsequent 
extensive theoretical works (e.g., Wambsganss 1990), many researchers focused on this 
interesting subject and presented many meaningful results. 

Roughly speaking, there are two different approaches to probe a structure of the 
central region of quasars. One is statistical approach by using long term monitoring data 
(e.g., 0stensen et al. 1996) that are expected to contain many microlens events. Recently, 
Wyithe et al. have performed thorough statistical study; they compare real and mock 
observational results, and constrain transverse velocities of lens, the mean mass of the lens, 
source size, and so on (Wyithe, Webster, & Turner 1999, 2000b, 2000c). Another is to 
focus on a single HME (High Magnification Event), on the basis of reasonable assumptions 
and/or some statistical features, and constrain a source size or structure (e.g., Wyithe et al. 
2000a, and references therein). 

In this letter, thus we take the latter approach and performed lightcurve fitting for 
the accurately and densely observed microlensing event that OGLE team have detected 
in Q2237+0305 lately (Wozniak et al. 2000a, 2000b), to put a limit on the microlensed 
source size. In section 2, we briefly present our method to fit an observed lightcurve and 
the results and discussions are shown in section 3. 



2. LIGHTCURVE FITTING 

The F-band monitoring data of OGLE team show that there is a dramatic brightening 
in image A and abrupt brightening and subsequent decay in image C (> 0.5mag), whereas 
flux changes have been less significant in images B and D (they have only ~ 20 % flux 
variation) compared with image A and C. This suggest that microlens events have occurred 
in images A and C independently. While the behavior of image A was complicated and 
could be caused by complicated (many) caustics, that of image C seems to be quite simple 
and can be understood in terms of a single-caustic induced HME. Therefore, applying 
approximate magnification formulae which are appropriate in the vicinity of a caustic, we 
try to fit the lightcurve of image C. Concrete formulae of magnification (/i) for a source 
position (x,y) are shown in Fluke & Webster (1999) for "fold caustic (hereafter, FC)" case 
and Zakharov (1995) for "cusp caustic (hereafter, CC)" case. 

To characterize properties of caustics, e.g., curvature of fold caustic, we should 
evaluate derivatives of Fermat potential (<f>) for gravitational lensing phenomena. In quasar 
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microlensing cases, Fermat potential is written as follows, 

= l{(v-x) 2 + ^-y) 2 -K c (r ] 2 + e)-l[(v 2 -^)cos2e + 2r ] ^sm2e}} 

- E s iog{Q[(^-^) 2 + (e-6) 2 ] 1/2 }, (i) 
i=i 

where, (rj, £) is an image position, 7 and 9 is shear and its direction, k c is convergence arising 
due to continuous (smooth) mass distribution, (77$, and in this equation represent 
the i-th lens position and its normalized mass, and iV lens is the total number of compact, 
stellar lens. Here, k c plus contributions of compact, stellar lens should be almost equal to 
total convergence (k) that is determined by macrolens model. All the length scales are 
normalized by the Einstein-ring radius (r E ). 

Unfortunately, we are not able to know complete spatial distribution and mass of 
compact, stellar lens objects, and it is impossible to characterize exact properties of caustics. 
Therefore, in this study, we assume that only one single lens object plays a significant 
role in the HME and contributions from other lens objects are negligible (Aq ens = 1), i.e., 
we only consider about so-called "Chang-Refsdal lens (Chang & Refsdal 1984)" situation 
(single lens object plus external convergence and shear), for simplicity, and characterize 
the properties. Applied values of convergence and shear is one of the best fit value on 
image C for macrolens model of Q2237+0305, k = 0.69 and 7 = 0.71 (Schmidt, Webster & 
Lewis, 1998). In some complicated case (e.g., lens objects are clustering), situation may 
be modified, or in the worst case, above assumption will break down. However, if this 
assumption holds (e.g., k c ~ k is the case), we are able to characterize the feature of FC 
and CC by only one parameter (#). 

In the case of the Chang-Refsdal lens, 2 different types of cusp caustic will be formed 
at 4 or 6 angular positions of the closed caustic curve (see Chang & Refsdal 1984), and 
the fold caustics appear as connecting these cusp caustics. Thus, we consider 2 cusp cases 
and 4 representative fold cases for an angle range of# = 0~7r/2 because of its symmetry. 
Moreover, we only consider about circular-shape source with top-hat brightness profile, i.e., 
neglect the effect of inclination angle, brightness profile. Therefore, magnification (M.) at 
any given source position (x, y) is obtained by the integration of /i on the circular disk with 
radius R. Assuming top-hat brightness profile is fairly simple treatment because the shape 
of microlensing lightcurve depends on source brightness profile (e.g., Yonehara et al. 1999). 
However, the resultant source size can be regarded as an effective (equivalent) source size 
and we use the top-hat brightness profile for convenience to numerical integrations (with 
1000 mesh number in this study). 



To include magnification caused by macrolensing and another caustics for microlensing, 



-4- 



here, we add constant magnification, M , plus gradual change of the magnification, M • t, 
to the total magnification. The latter mimics magnification changes for the ensemble of 
other microlens and this gradual change make the fits better (Wyithe and Turner, private 
communication). Furthermore, intrinsic variability of this quasar may be long duration 
with small amplitude (Wyithe et al. 2000d), and this term may also mimic the intrinsic 
variabilities of quasar. 

Additionally, we should evaluate the apparent magnitude of quasar without any 
microlensing and macrolensing, i.e., the intrinsic magnitude, tuq. It is quite difficult because 
quasars have intrinsic variabilities and also the observed magnitude of Q2237+0305 is 
affected by small amplitude and/or long timescale microlensing. Fortunately, from the 
monitoring data by 0stensen et al.(1996) and the OGLE team, the observed magnitude 
of image C is roughly constant at ~ 18.6 mag long before the recent HME. We thus take 
this as the magnitude of this quasar without any microlensing. Moreover, by using the 
previously applied value of k and 7, we can evaluate the macrolens magnification of image 
C, and the intrinsic magnitude is estimated to be mo = 18.6 — 2.5 log [(1 — k) 2 — 7 2 ] ~ 19.6. 

Finally, assuming the source trajectory is straight and determine source velocity on the 
source plane, v = (v x ,v y ), time when the source crosses a caustic (FC case) or x-axis (CC 
case), To, the expected microlensing lightcurve for any given parameter is obtained from 
m(t) = m — 2.5 log M + M ■ t + A4(t) (for CC case and fold caustic gazing case, impact 
parameter, d is also required). 

To obtain the best fit lightcurve and its parameter, we minimize x 2 value between the 
observed lightcurve, m bs(£)j and the mock lightcurve for a given parameter, m(t), for each 
caustic case (FC and CC) by using one of downhill simplex method, so-called "AMOEBA" 
routine (Press et al. 1986). Nonetheless, caustic is a kind of singularities, and fitting 
methods including singularities do not work so well. For this reason, we subdivide each 
caustic case into all considerable path cases (depicted in figure P and perform lightcurve 
fitting at every possible case. After the best fit parameters are obtained, we compare reduce 
X 2 for all the possible cases for each of FC and CC, and determine the best fit (smallest x 2 ) 
parameters for each FC and CC case. 

In this lightcurve fitting, we only took into account data points around the peak of 
image C (JD - 2450000 = 1289.905 ~ 1529.531, total number of used data points is 83). 
Of course, we can also taken into account data points before and after the peak, but in 
those epochs, distance between the source and caustic could be larger than that in the peak 
region. Thus, it is not clear whether the approximation for magnification is reasonable or 
not before and after the peak. And so, we have restricted data points only around the peak. 
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3. RESULTS AND DISCUSSIONS 

Our resultant, best fit parameters for all the cases that we have considered are 
summarized in table [I]. Some resultant path of source, the best fit lightcurve, source size 
dependence of reduced x 2 are also shown in figures || and || (degree of freedom is 83 — 5 = 78 
for FC, 83 — 6 = 77 for CC). In the case of a infinitely-small size source, expected lightcurves 
for microlens event are quite different from case to case, i.e., a lightcurve for FC case and 
that of CC case is clearly different, and also, difference between that of fold-1 and fold-2 
(figure [1]) is evident and so on. But, as you can easily see in these table and figures, the 
fitting for all cases works well (the best-fit reduced \ 2 is ~ 1)- This fact owing to the 
finite-size source effect, and we can manage to reproduce the observed feature at every 
considerable case. 

We also performed a Monte-Carlo simulation for every case to estimate confidence 
region of best fit parameter by using following procedures; (1) Supposing that the best fit 
parameter is real parameter, we calculate an ideal lightcurve without any errors, (2) Add 
random errors with the magnitude corresponding to the observational error dispersion, and 
sample this lightcurve at the times corresponding to the actually observed times, we obtain 
a mock lightcurve. (3) By using this mock lightcurve, we perform a lightcurve fitting again 
for all considerable cases (at FC and CC case indicated before, see also fig [I]) and obtain a 
set of the best fit parameters for the mock lightcurve. 

Iterating procedure (2) and (3) for 100 times in this study, summarizing the best fit 
values for mock lightcurves, and we can evaluate a confidence region. To evaluate the 
90% confidence region from Monte-Carlo results, we calculate total \ 2 between lightcurve 
which is actually observed and that is obtained from the parameters of Monte-Carlo result. 
Subsequently, we pick up 90% parameter sets which have smaller total \ 2 - From these 
selected parameter sets, finally, we can obtain maximum and minimum values of parameters 
and we define ranges between these maximum and minimum as the 90% confidence region. 

In figures and |3|, we also presented a histogram for the x 2 differences between "mock" 
lightcurves and the best fit lightcurve, and ideal \ 2 distribution curves for corresponding 
degrees of freedom (5 for FC, 6 for CC). These two exhibit similar distributions and our 
confidence region estimate seems to be reasonable. 

For every case, if the source size is larger than the best fit value, expected magnification 
will be suppressed, lightcurve will become shallow, and goodness of fit is reduced. On the 
other hands, if the source is smaller than the best fit value, expected magnification will 
be enhanced, lightcurve will become sharp, and goodness of fit reduced, too. These are 
qualitative reason why the source size is limited in somewhat small range. 
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Considering all our fitting result, at least, we can say that the source size of Q2237+0305 
should be smaller than ~ 0.98 Einstein-ring radius (more than 90% confidence level). This 
upper limit is given in the case of FC (fold 3 in table [I]), while another case suggests a much 
smaller source size. Thus, this limit is a fairly conservative upper limit for the source. Since 
the Einstein-ring radius of quasar microlensing is typically £ 1 /xas, our result indicates the 
existence of a sub-/xas source in quasar ! 

To obtain the actual size, we have to calculate r E for relevant parameters. If we assume 
1.0M© as a mass of lens object (Mi ens ) and the Hubble constant Hq = 67km s" 1 Mpc -1 
(Kundic et al. 1997), r E will be corresponds to 10 17 cm. This value strongly depends on H 
(rg oc H 1 ^ 2 ) and the lens mass (rg oc rather than the cosmological parameters in 

this case (~ 10% uncertainty covers roughly all the reasonable range). And finally, we get 
10 17 cm ~ 7 x 10 3 AU as the resultant upper limit to the source size. This value is consistent 
with the result of the statistical research performed by Wyithe et al. (2000c). Moreover, 
for the ~ O.1M lens object case suggested by Wyithe et al. (2000b) as a mean lens mass, 
the size will be reduced by some factor and become ~ 2 x 10 3 AU !. Alternatively, the 
lens object may be a stellar object, and so, there is an upper mass limit to exist stably. 
Even if we adopt this upper limit ~ 1OOM for the lens mass, the size will be ~ 0.3pc, at 
most. Therefore, our result strongly supports the existence of an accretion disk in a quasar 
and that the accretion disk smaller than this size is a fairly dominant source of radiation 
from the quasar, at least in the F-band at observer frame. Additionally, resultant effective 
transverse velocity on the source plane is £ 10 5 (km s _1 ), and is also consistent to the value 
presented by Wyithe, Webster & Turner, 1999. 

On the other hands, if we assume that the accretion disk is a type described by 
the standard accretion disk model (Shakura & Sunyaev 1973), we can also estimate 
the effective source size from its luminosity. The magnitude of this quasar in the 
absence of a macrolens effect is easily converted into flux, f v (Wozniak et al. 2000a). 
If we denote the absorption as Ay mag and adopt the luminosity distance (g?l) to the 
quasar, luminosity of the quasar at this observed waveband (L) will estimated to be 
L ~ is(fv ■ 10°- 4Av )47tc?l ~ 3.7 x io 42+a4A v er g s -i_ Furthermore, radiation process of the 
standard accretion disk is blackbody radiation, and the effective temperature (T efT ) of the 
accretion disk at this waveband corresponds to ~ 2.5 x 10 4 K. On the other hands, we can 
relate the effective temperature, luminosity and the radius through the central black hole 
mass and accretion rate, 

Mi)" Vvi/2 - »> 

where, a is the Thomson-scattering cross section. Consequently, we are able to estimate the 
effective source size being r ~ 2.0 x iq 14 +°- 2A v cm There is an uncertainty in Ay, but this 
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is consistent with our results and strongly indicates the existence of an accretion disk. We 
should pay attention that we do not insist the existing source should be a standard-type 
accretion disk (other accretion disk models may be also consistent with our results). Rauch 
& Blandford (1991) reported that the accretion disk in Q2237+0305 is non-thermal or 
optically thin. But, there are some ambiguities in this work (e.g., absorption, lens mass), 
and it is quite difficult to support or oppose the report from our results. To specify disk 
models, we have to do more extensive monitorings and/or analyze multi-band microlens 
lightcurves to compare the resultant source sizes. 

The performed fitting procedure do work, but strictly speaking, all our best-fit reduced 
X 2 are somewhat larger than 1. This means the goodness of fit of our results is not extremely 
good and probably, there may be some systematic errors that we do not take into account. 
In this work, we neglect the detailed feature about source, magnification patterns and 
intrinsic variabilities of the quasar. However, if such effects are really existing, the shape of 
the lightcurve will be systematically altered and the best-fit reduced x 2 ma y be increased 
by these effects. There are difficulties to be take into account all above possibilities in our 
procedure, but it will be done in future. 

In future work, we should further develop quasar microlens technique in two statistical 
ways. One is to study statistical properties of magnification near FC and CC. Although 
the statistical features have already been studied by many researchers (e.g., Wambsganss 
& Kundic 1995), the effect, such as, lens object clustering that affect properties of 
magnification in the vicinity of caustic is not well understood. The other is to continue 
monitoring this kind of quasars, sample similar microlensing events more and more. Such 
analysis may be able to reduce ambiguities arising due to unknown lens mass, and/or 
different features of caustic networks and so on. 
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valuable suggestions and comments. The author also acknowledge the OGLE team for 
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Table 1: Best fit parameters for several possible HME (see figure p]) and their reduced \ 2 (x 2 )- 
All the length scales and time scale is normalized by and one day. The unit of velocity is 
the Einstein-ring radius divided by a day. For T , T = correspond to JD — 245nnnn = n. 
Upper and lower value denoted beside the best fit parameters show a 9n% confidence level 
which calculated from Monte-Carlo simulation (see the text). 
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Fig. 1. — Schematic view of possible HME. All the six cases to reproduce the observed, fairly 
symmetric HME event are shown: 2 are fold caustic cases and 4 are cusp caustic. 
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goodness of fit results of Monte-Carlo 




Fig. 2. — Fitting results in a case of fold caustic (fold-3). In upper left panel, the best fit 
source path (dashed line) and relative to caustic (bold line) is depicted. Upper right panel is 
observed flux (error bar) and the best fit lightcurve (bold line). Lower left panel is total x 2 
value respect to source size (bold line), and result of Monte-Carlo simulation (histogram). 
Lower right panel distribution of x 2 with 5 degree of freedom (bold line) and distribution 
of x 2 difference between mock lightcurve and the best fit lightcurve (histogram). All the 
length scale is te- Kinks in source size dependence of total x 2 (lower left panel) are caused 
by changes from a best-fit subdivided case to another subdivided case, e.g., at source size 
~ 0.6, the best-fit, subdivided case changes from FC1 to FC2. 
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Fig. 3. — Same as figure |2| but cusp caustic case (cusp-1). Degree of freedom of lower right 
panel is 6 in this case. Kinks in source size dependence of total \ 2 value (lower left panel) 
appeared with the same reason as in figure § 



